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Abstract We analyse the possibility of distinguishing Dirac and Majo- 
rana neutrinos in future neutrino factory experiments in which neutrinos are 
produced in muon decay when, in addition to a vector type as in the SM, 
there are also scalar interactions. We check this possibility in an experi- 
ment with a near detector, where the observed neutrinos do not oscillate, 
and in a far detector, after the neutrino oscillations. Neglecting higher-order 
corrections, even neutrino observation in the near detector does not give a 
chance to differentiate their character. However, this possibility appears in 
the leading-order after the neutrino oscillations observed in far detector. 

1 Introduction 

In the Standard Model (SM) with only V-A coupling between neutrinos 
and charged leptons, massless Dirac and Majorana neutrinos are not distin- 
guishable. In the case of massive neutrinos, whose masses are many orders 
of magnitude smaller than the masses of charged leptons, the distinction be- 
tween these two types of neutrino is formally possible but very difficult [T]. 
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Indeed, the experiments so far have not given us any guidance on the question 
of whether neutrinos are Dirac or Majorana particles. Probably neutrino-less 
double beta decay process, if observed, will be the first experiment that could 
indicate whether the neutrinos are Majorana in nature. When neutrino inter- 
actions have also non-standard contributions, there are other possibilities for 
distinguishing between Dirac and Majorana neutrinos. We will discuss one 
such possibility in this Letter, related to muon decay. However, the situation 
is not so simple. In paper [2], muon decay was analysed assuming the most 
general local, derivative-free, Lorentz invariant and lepton nonconserving in- 
teractions. Assuming that neutrinos are not observed and their masses can 
be neglected, it was shown that up to the leading order in the parameters 
beyond the SM, it is impossible to distinguish between Dirac and Majorana 
neutrinos. In this Letter we will go further and analyse the problem of deter- 
mining whether neutrinos are Dirac or Majorana particles when muon decay 
is used as a source of neutrinos in oscillation experiments (future neutrino 
factories). We consider neutrino detection in the near (no oscillation) and 
far (after oscillations) detectors. 

In the next section we will consider Dirac and Majorana neutrinos oscil- 
lation in the case of their non-standard interaction. Section 3 contains our 
conclusions. 

2 Oscillation of Dirac and Majorana neutri- 
nos 

We assume that neutrinos are produced in muon decay as is the case for 
future neutrino factories [3j. In the SM, muon decay can be described by an 
effective Lagrangian 

Li = -2V2Gf [gl {ei'^PLvi) (I^7«^l/^)] + h.c, (1) 

with related to the elements of the neutrino mixing matrix gj^ = UeiU*^. 
However, if we go beyond the SM, different Lorentz structures can appear. 
Present experiments give bounds on the effective couplings for different types 
of interactions, i.e., scalar vector and tensor, see, e.g., |ll[5]. A detailed analy- 
sis of these constraints shows that from among all possible types of couplings, 
the scalar one can give the most significant contributions. Therefore we con- 
fine ourselves to this additional coupling and we take our Lagrangian to be 
of the form 

Ci = -2V2Gf [4 (ePRUi) (ujPLf^) + 4 (ct^Plz/.) O^laPLft)] + h.c. (2) 
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Many SM extensions predict additional interactions that contribute to muon 
decay and generate the effective Lagrangian presented above (see, e.g., [6l 
[3, IH]). Another reason for choosing this type of Lagrangian is that it is 
the only type of interaction where for Majorana neutrinos with negligible 
mass, scalar and SM vector amplitudes interfere even in the limit of the 
electron mass going to zercRgiving the opportunity for leading order effects 
[9]. The matrix in Eq.pl) can be chosen freely in agreement with the 
present bound, while for the vector interaction we want to have it the same 
value as in the SM so we choose it such that it is proportional to its SM 
value i.e. {g^)'^g^ = Wg^W^Pn, where is a projection operator in neutrino 
flavour space, subscript /i means that it projects on muon neutrino type 
and \\g^\\ = \/Tr[{g^y g^] is the Hilbert-Schmidt norm. In mass base the 
projection operator is given by {Pfj.)ij = Uf^iU*j. This is equivalent to 
assuming that in our case 

4 = <7ef/«f/;„ (3) 

where gc is a constant, expected to be close to 1. A more complicated vector 
interaction than (3), affecting flavour symmetry conservation, has been con- 
sidered in the literature (see, e.g., [H HO]), but without any impact on our 
overall conclusions. 

It was noticed [2] that in muon decay, in the case when the neutrinos are 
not directly observed and have negligible masses, Dirac and Majorana neu- 
trinos are impossible to distinguish. We will show that this conclusion holds 
also for the case of a near detector but it is no longer true when neutrinos 
oscillate. We will use the density matrix formalism for describing neutrino 
states, since this properly takes into account the interference between dif- 
ferent amphtudes. The usual treatment based on pure states [H] cannot be 
applied here [121 1131 E] • The density matrix formalism which we adopt will 
enable us to describe neutrino states as proper mixed states in the quantum 
mechanical sense, giving information about the neutrino spectrum produced 
in muon decay. We briefly describe this formalism (more detailed information 
can be found in [El [131 H US]). 

We denote by A{X,i, E,lj) the amplitude describing neutrino production 
in muon decay in the mass state i, helicity A, and energy E; u characterises 
all discrete and continuous degrees of freedom of the other particles that 
are produced with the neutrinos. Then the matrix elements of the density 

^In subsequent calculations we will neglect the electron mass. Non-vanishing electron 
mass would introduce another possibility for interference but this effect will be suppressed 
by powers of the ratio of electron mass to muon mass ^ which is below the expected 
accuracy of future neutrino factories. 
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operator ar^ 



{Q{E))x,i,x'i' = ) [V 

where is a normalisation constant chosen so that 



J Ern.i.n. \ 



and our density matrix contains information about the state of the neutrino 
as well as its energy spectrum. In our calculations we assume that the initial 
muons are not polarised, so we sum also over the muon polarisation states. 
Straightforward but tedious calculations lead to the following result. The 
state of a Dirac neutrino produced in muon decay, neglecting the mass of the 
electron, is 



_ 4x\2{g^yg^{3 - 2x)P_ + 3(/)t/(l - x)P+) 
Tr[{gSygS]+4Tr[{gyygy] ' 

while in the case of the Dirac antineutrino, we have 

_ 2x\24{g''f{g^y{l - x)P+ + (/)^(/)*(3 - 2x)P. 
^ Tr[{gS)^gS]+ATr[{gy)^gy] 

and for a Majorana neutrino we obtain 

irtir(3 - 2x)x^P^ + 6K*K^{1 - x)x^P+ 



Tr[K^K] 



(5) 



(6) 



(7) 



where the K matrix is given hy K = g^ + \{g^)^, x = ^ {E is the en- 
ergy of the neutrino and M is the muon mass), and P^ (respectively, P_) 
is the positive (respectively, negative) helicity projection operator. The nor- 
malisation is such that dxTr[g\ = 1. Since we assume that the initial 
muons are not polarised, the neutrino density matrices do not depend on 
the neutrinos' direction of flight. These formulas (Eqs.(5-7)) are given in the 
muon rest frame and the lack of manifest Lorentz invariance in our formal- 
ism might raise some doubts. Fortunately, for very small neutrino masses, 
Lorentz transformations do not affect the mass-spin structure of the neutrino 
density matrix (see, e.g., [T2]). 



^In principle, in our treatment, the density matrix elements should also be labelled by 
a pair of continuous indices E,E', however using energy-momentum conservation it can 
be shown that density matrix is diagonal with respect to energy. So we will omit the delta 
function 6{E — E') and treat the density matrix as a scalar function of energy. 
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For simplicity, in what follows we will assume that in the detection pro- 
cess we measure only the neutrinos with one specific helicity, e.g., left-handed. 
This assumption does not restrict the generality of our considerations, since 
any detection process can be represented by an incoherent mixture of pro- 
jection operator on both neutrino helicity state^ In the case when only one 
helicity state is detected, the most general operator describing the detection 
process is proportional to the projection operator on that one specific helic- 
ity. Using the properties of projection operators, we conclude that all the 
energy dependence of the density matrix can be factored out in such a case. 
This energy dependence affects the flux of the neutrinos. Since we are inter- 
ested in the difference between Dirac and Majorana neutrinos and not the 
total number of predicted events, this part is not important for our study. 
In our case, the factor (3 — 2a;) is absorbed into the flux of the neutrinos 
both for the Dirac and the Majorana case. Another important point is that 
the coupling matrices and must be normalised in such a way that Eq. 
(|2]) leads to a value of the decay rate which is in agreement with precision 
measurements of this observable. This means that for Dirac neutrinos we 
must include condition [3] 

lrr[(/)t/] + Tr[(^7^)t^n = l- (8) 
For the analysis which follows, we can rewrite this condition as 

ii/ii = /^^^w ^ 1 - ^ (9) 

In the Majorana case, the condition is different [2]. Due to the interference 
between the scalar and vector amplitudes, 

Tr\K'^K\ = 1. (10) 

This can be solved for the norm of the matrix 

ll/ll = i/2(v/4H¥W^-«II/II) (11) 

- i-^ + ^(«^-i)ii/ir+o(ii/ir), (12) 

with a = Re{Tr[gff{g]^)*]) , g^ being a matrix normalised so that its norm is 
unity, and similarly for gf^. We observe that for a = 0, i.e., no interference 

•^This basically means that we can measure only left- and right-handed neutrinos and 
not a linear combination of them. 
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between the scalar and vector contributions, this condition is the same as in 
the Dirac case. 

Adding this up, we can now consider for the Dirac case a neutrino state 
described by the matrix 



QD = {g'^Yg'', (13) 



and for the Majorana case. 



0M = K^K, (14) 



with and chosen such that Eqs.(|8|and lOl) are fulfilled. The oscillation 



is described as usual by the unitary transformation 

q,{L) = U[L\q,U[L\\ x = D,M, (15) 

where U[L] sets the neutrino propagation in vacuum or in matter, depending 
on the experimental conditions. In general, in the case of vector right-handed 
interactions, the propagation in matter can also lead to the possibility of dis- 
tinguishing Dirac and Majorana neutrinos [16]. However we are considering 
in this Letter only a scalar interaction, so we assume that the propagation 
Hamiltonian is the same for the Dirac and Majorana neutrinos. In numeri- 
cal studies, we assume that the neutrinos propagate in vacuum, however the 
generalisation to matter is straightforward and changes only the quantitative 
results without changing the general pattern. 

Let us then calculate the number of observed neutrinos. With our as- 
sumptions neither the detection cross-section for a typical detection process 
p!7] nor the flux depends on the nature of the neutrino, so we can only 
consider the oscillation probability 

(L),. =Tr[&[L]P«], x = D,M. (16) 

Here, Pa is a projection operator in neutrino flavour space, which projects 
to flavour a direction, {Pa)i,j = UaiU*j. 

We neglect the NP contribution to the detection process because it cannot 
depend on the nature of the neutrino when no vector right-handed interac- 
tions are present. 

Calculating the probability factor (16) for L = using conditions (|9|, we 



obtain for the Dirac case 

-,5'||2 



P,MO)n = 1-^, (17) 
Pw(0)d = 0. (18) 
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For the Majorana case, using (11), we obtain 



II 5112 

P,^,{0)m = 1-^^(1-U. (19) 

P^^e,r(0)M = ^^ee,r. (20) 

The = Tr[{gff)* (gff)'^ Pa] depends only on the flavour structure of the 
matrix and can take any value between and 1. The difference between 
the Dirac and Majorana probabilities for L = is of second order in the 
strength of the scalar interaction: therefore one expects that this difference 
is negligible, so we can conclude that in practice there is no way to distinguish 
Dirac and Majorana neutrinos at the near detector. It is usual, with non- 
standard interactions, to take into account only effects that are linear in the 
non-standard parameters and neglect the higher order corrections |llj . 

The key point in this calculation was the normalisation condition (Eq.(|To|)), 
from which the linear term for the Majorana neutrino disappears, as then 
the dependence on the New Physics strength is the same for both types of 
neutrinos. This no longer occurs when we allow the neutrino to oscillate. Let 
us only keep the terms linear in the norm of the scalar interaction. Then a 
simple calculation for Dirac neutrinos gives (/3 = e,yU, r) 

P,ML)d = P,ML)sM + 0{\\g'\\^), (21) 

with Pfj,^j3{L)sAi the oscillation probability calculated assuming only the SM 



contribution. For the Majorana case, from Eqs. (11), (14), (15), and (16), 
we obtain 

p,ML)m = \\g''\\'P,ML)sM + Wg^^W ■ \\g'\ML) + 0(||/in (22) 
= P,ML)sM - \\g'\\{aP,^^{L)sM - ap{L)) + 0(||/|n, 

where afs{L) = Re{Tr[U[L]{g^y{g]^fU[L]^P,3])- While obviously the prob- 
ability factor for Dirac neutrinos still depends quadratically on the strength 
of the scalar interactions, for the Majorana neutrino the dependence is now 
linear. It must be kept in mind that the parameters have different 

meanings for Majorana and Dirac neutrinos. This however does not prevent 
us from concluding that Dirac and Majorana neutrinos can have different os- 
cillation rates. We have also checked this numerically, by generating random 
g^ matrices and calculating the probability factor. The results are given in 
Fig. 1, which shows the probability P^^e{L) as a function of L/ E[km/GeV] 
for \\g^\\ < 0.1 for both Dirac and Majorana neutrinos. For the Dirac case, 
the effects of the neutrino scalar interactions are practically indistinguishable 
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Figure 1: Pf^-^e{L) as a function oi L/E for Dirac and Majorana neutrinos. 
The matrices were randomly generated satisfying < 0.1. 

from the SM contribution. However for Majorana neutrinos the effects are 
larger. 

Fig. [2] shows P^^^(L) for L/E = lOOOkm/GeV as a function of 
Current experimental limits give \\g^\\ < 0.55 [1]. From this plot it is clear 
that observation of the deficit of muon neutrinos on the level outside the 
allowed range for Dirac neutrinos would indicate that neutrinos are Majorana 
particles. However, if we observe a number of neutrinos that is in agreement 
with the current limits for Dirac neutrinos then we cannot say anything about 
whether neutrinos are Dirac or Majorana particles. 



3 Conclusions 



We calculated the oscillation probability for Dirac and Majorana neutri- 
nos in the presence of additional scalar interactions. The dependence on the 
strength of the new interaction is different for Dirac and Majorana neutri- 
nos. The appearance of a linear term in the expansion for the Majorana 
neutrinos can in principle help to determine whether the neutrino is a Dirac 
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Figure 2: P^^^(L) for ^ = 500j^ as a function of ||5''^||. 

or a Majorana particle. Unfortunately, this analysis is model-dependent and 
if the deviation were within the allowed range for Dirac neutrinos, then we 
cannot say anything about the nature of the neutrino. When no oscillations 
are present, i.e., for L = 0, in the first order for the transition probabilities 
the dependence on the strength of the scalar interaction is the same for Dirac 
and Majorana neutrinos. 

Even though these kind of predictions are model dependent and not al- 
ways conclusive, they can hint at some New Physics and at the determination 
of the nature of the neutrino as a byproduct of planned experiments, such as 
a neutrino factory. 

References 

[1] L. F. Li, F. Wilczek, Phys. Rev. D25 (1982) 143; B. Kayser, R. E. 
Schrock, Phys. Lett. B112 (1982) 137; B. Kayser, Phys. Rev. D26 (1982) 
1662. 

[2] P. Langacker and D. London, Phys. Rev. D 39 (1989) 266. 



9 



[3] S. Geer, Phys. Rev. D57 (1998) 6989, ibid. 59 (1999) 039903E; A. De 
Rujula, M. B. Gavela, and R Hernandez, Nucl. Phys. B 547 (1999) 21. 

[4] J. Beringer et al. (Particle Data Group), Phys. Rev. D 86 (2012) 010001. 

[5] F. del Aguila, J. de Bias, R. Szafron, J. Wudka, and M. Zralek, Phys. 
Lett. B 683 (2010) 282. 

[6] P. Duka, J. Gluza, and M. Zralek, Annals Phys. 280 (2000) 336. 

[7] W. Fetscher, H. J. Gerber, and K. F. Johnson, 1986, Phys. Lett. B 
173, 102; Y. Kuno and Y. Okada, Rev.Mod.Phys. 73 (2001) 151-202; 
S. Antusch, J. P. Baumann, and E. Fernandez- Martinez, Nucl. Phys. B 
810 (2009) 369. 

[8] M. Medina and P. C. de Holanda, Adv. High Energy Phys. 2012 (2012) 
763829. 

[9] R. Szafron and M. Zralek, Acta Phys. Polon. B 40 (2009) 3041. 

[10] P. Coloma, A. Donini, J. Lopez-Pavon, and H. Minakata, JHEP 1108 
(2011) 036. 

[11] Y. Grossman, Phys. Lett. B 359 (1995) 141. 

[12] M. Ochman, R. Szafron, and M. Zralek, J. Phys. G G 35 (2008) 065003. 

[13] R. Szafron and M. Zralek, Prog. Part. Nucl. Phys. 64 (2010) 210. 

[14] R. Szafron and M. Zralek, Acta Phys. Polon. B 42 (2011) 2501. 

[15] M. Ochman, R. Szafron, and M. Zralek, Nucl. Phys. Proc. Suppl. 217 
(2011) 347. 

[16] F. del Aguila, J. Syska, and M. Zralek, Phys. Rev. D 76 (2007) 013007. 

[17] A. Ccrvcra, F. Dydak, and J. Gomez Cadenas, Nucl. Instrum. Meth. A 
451 (2000) 123. 



10 



